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Particle sorting by a structured microfluidic ratchet device with tunable selectivity: 

Theory and Experiment 

Lukas Bogunovic, 1 Ralf Eichhorn, 2 Jan Regtmeier, 1 Dario Anselmetti, 1 and Peter Reimann 1 '!^ 

1 Universitdt Bielefeld, Fakultat fur Physik, 33615 Bielefeld, Germany 

B NORDITA, Roslagstullsbacken 23, 106 91 Stockholm, Sweden 

(Dated: March 1, 2013) 

We theoretically predict and experimentally demonstrate that several different particle species 
can be separated from each other by means of a ratchet device, consisting of periodically arranged 
triangular (ratchet) shaped obstacles. We propose an explicit algorithm for suitably tailoring the 
externally applied, time-dependent voltage protocol so that one or several, arbitrarily selected parti- 
cle species are forced to migrate oppositely to all the remaining species. As an example we present 
numerical simulations for a mixture of five species, labelled according to their increasing size, so 
that species 2 and 4 simultaneously move in one direction and species 1, 3, and 5 in the other. The 
selection of species to be separated from the others can be changed at any time by simply adapting 
the voltage protocol. This general theoretical concept to utilize one device for many different sorting 
tasks is experimentally confirmed for a mixture of three colloidal particle species. 
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I. INTRODUCTION 

Ratchet effects refer to directed transport under non- 
equilibrium conditions in periodic systems with broken 
spatial symmetry, and are currently attracting consider- 
able attention in a wide variety of different fields [l|-(j]- 
So far, the majority of pertinent works has been devoted 
to basic theoretical concepts, fueled - among others - 
by the exciting perspective to exploit ratchet effects for 
particle sorting purposes at the nano- and micro-scale. 
Concerning experimental ratchet devices, directed trans- 
port per se has been demonstrated by now in quite a 
number of systems as well, see refs. |5Hl7| for just a 
few examples. On the other hand, the actual sorting 
of different particle species has been achieved only in a 
few experimental systems [18-29]. All of them have been 
conducted in some microstructures with broken spatial 
symmetry, while the indispensable non-equilibrium con- 
ditions have been implemented mostly by utilizing elec- 
trokinetic effects [18l - [27| , but also via magnetic fields [28( 
or intracellular morphological changes [23 ] . 

Along these lines of experimentally realizing ratchet 
effects by exploiting electrokinetics [30, [31| , our present 
work addresses the natural next question: Given a mix- 
ture of different particle species, is it possible to separate 
any arbitrarily selected subset of species from all the oth- 
ers with the help of one and the same microstructured 
ratchet device? After setting the general framework we 
have in mind, we will outline the basic theoretical ideas, 
verify them quantitatively by numerical simulations, and 
finally provide the experimental proof-of-principle for a 
colloidal suspension containing three different species of 
polystyrene beads. 

The general experimental framework is schematically 
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FIG. 1: Experimental setup, (a): Schematic top view (x- 
y-plane) of the entire device (not to scale) with reservoirs 
and electrodes at its ends and a central, microstructured 
part (grey), (b) Optical micrograph image of the microstruc- 
tured area. The entire microstructured region extends over 
10000 (imx 300 fim in the x-y-plane and 10 /nn in height (z- 
direction). (c) Magnification, exhibiting periodically arranged 
obstacles (posts) of ratchet-shaped cross section (see Fig. 2 for 
the dimensions of the obstacles) with a period of 20 /im in x di- 
rection and 14 /jm in y direction. Also visible are microbeads 
of three different diameters (2.9 /im, 1.9 /jm, and 1.1 ^m). The 
largest and smallest beads can be easily distinguished by their 
size. The medium sized beads are fluorescently labelled and 
thus appear as bright spots. 
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illustrated in Fig. 1 (see also Appendix A.l). The central, 
topographically structured area consists of periodically 
arranged, triangular (ratchet) shaped obstacles. The en- 
tire device is filled with a buffer solution, containing par- 
ticles in sufficiently high dilution, so that their mutual 
interaction is negligible. It is actuated by an external 
voltage like indicated in Fig. la. 
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FIG. 2: Force fields, (a) Electrical field E\[f) around one 
idealized obstacle within a region of size 20 /jm x 14 ^im, cor- 
responding to one basic "unit cell" of the periodic microstruc- 
ture, see also Fig. lc and Eq. ([2]). The shape of the obstacle is 
chosen to closely approximate the real experimental shape of 
the triangular posts as formed by the soft lithography process 
(see Appendix A.l) and as characterized by REM imaging. 
The obstacle extends 8.5 /im from baseline to tip (a; direc- 
tion) and 8.5 /im from side vertex to side vertex (y direction), 
(b) Dielectrophoretic force field aVfii(r)/2, under the as- 
sumption that q < 0, see also Eq. ((5]) andAppendix A. 5. 
Red dots indicate particle traps (minima of the potential ([5])). 
Both fields inherit the spatial periodicity of the microstruc- 
ture (Fig. lb) and are depicted here in arbitrary units. 



II. THEORETICAL CONCEPTS 

A. Electrokinetic effects 

Upon applying an electrical DC- and AC-voltage U(t) 
of the general form 



e, its viscosity 77, and the particle's zeta potential C by 
the well-known Hclmholtz-Smoluchowski relation 



U(t) = Uuc + U AC sin(wi) 



(1) 



to the electrodes in Fig. la, an electrical field E(r, t) is 
created at position f = (a;, y, z) and time t within the 
device. For all experimentally realistic DC- and AC- 
voltages in (p}, the electrostatic approximation is valid 
[321 ], so that the field induced by U{t) is 



E{f,t) = U(t)E 1 (r)/U 1 , 



(2) 



where E\{r) denotes the electrical field induced by a 
unit DC-voltage U\ :— IV. To determine Ei(r) within 
the structured region (Fig. lb), we numerically solved 
Laplace's equation with Neumann boundary conditions, 
thus approximately treating the microstructure as a per- 
fect isolator [32|, [33| . The resulting electrical field E\ (r) 
around an experimentally realistic obstacle is illustrated 
in Fig. 2a. Outside the structured region, an approxima- 
tive treatment along the lines of [34j was adopted to es- 
timate the voltage drop over the structured region, given 
the voltage U\ := IV is applied to the device electrodes 
(cf. Fig. 1). 

Due to the electric field E(f,t), a particle in the mi- 
crostructure experiences two different electrokinetic ef- 
fects [3l|, |35[: First, it is driven by electrophoresis. The 
particle velocity r due to purely electrophoretic forces 
(i.e. neglecting other electrokinetic effects, thermal noise 
etc.) can be expressed in terms of the fluid's permittivity 



r = E(f,t)eC/v, 



(3) 



provided the Debye screening length of the counter-ions 
is much smaller than the particle diameter [35|, [36] . The 
latter condition can be safely taken for granted in view 
of a Debye length of ~ 10 nm for our experimental 10 
mM phosphate buffer [35[ . The Helmholtz-Smoluchowski 
relation ^ applies to particles of any shape [36| and 
holds also if the external electric field varies over dis- 
tances comparable to the particle size [37j. We further 
note that it can formally be rewritten as a force balance 
•yr = q*E(f) between the viscous (Stokes) friction force 
7 r with friction coefficient 7 and an effective electrostatic 
force q*E(r) with an "effective charge" q* defined as 

q*:=l<h (4) 

This quantity q* is a purely formal abbreviation which 
has the dimension of a charge but has very little to do 
with any kind of actual physical charge. 

Second, the electrical field E(r, t) induces a dipole mo- 
ment p in the particle and the ion cloud around it, which 
is characterized by the relevant (effective) polarizability 
a. In general, a is a complex, frequency dependent quan- 
tity, and the dipole moment p may be composed of sev- 
eral contributions according to the frequency components 
contained in the electrical field E(r, t). Under our present 
experimental conditions, however, the quasi-static ap- 
proximation with a referring to the zero frequency limit 
can be invoked, since polarization processes are much 
faster than time- variations of the electric field [38[ • Like- 
wise, conductive and dielectric losses can be safely ne- 
glected, so that a is given as a real-valued quantity. The 
dipole moment p is thus of the form a.E{f,t), and the 
standard dipole-field-coupling force (p ■ \J)E derives from 
a dielectrophoretic potential 



W(r,t) = -aE\r,t)/2 



(5) 



Figure 2b depicts the corresponding dielectrophoretic 
force field around an experimentally realistic obstacle. 
Similar as q* , also the polarizability a is an effective 
quantity depending on particle and solution properties. 
A common approximation for the present situation of a 
colloidal bead (diameter d) suspended in an ionic buffer 
solution is 



a = -~2ire(d/2) 3 



(6) 



obtained from the well-known basic formula for the po- 
larization of a non-conducting dielectric sphere in a con- 
ducting dielectric medium with permittivity e at low fre- 
quencies of the electric field [38| . 



B. Particle dynamics 

In addition to the two electrokinetic effects from the 
previous Section, other relevant forces acting on the par- 



tide are thermal fluctuations and interactions with the 
microstructure, while inertia effects are so tiny that they 
do not play any noticeable role [3a, |39| • We model the 
particle dynamics by the standard ovcrdamped Langevin 
equation [l-|I [UlI [H, [H [H [H 

7 r* = q*E(f, t) + aVE 2 (f , t)/2 + F{r) + ({t) . (7) 

On the right hand side, the first and second terms are the 
electrophoretic and dielectrophoretic forces as detailed 
above. The third term F(r) accounts for the repulsion 
by the microstructure, idealized as hard wall forces due 
to the periodic obstacles from Fig. 2 and the channel 
walls. The last term £(£) describes thermal fluctuations. 
Its components £<t(£), <j £ {x, y, z}, are modeled as usual 
[l|-(j] by independent Gaussian noises of zero average and 
correlation (£ CT (t) £ CT (s)) = 2jkT5(t — s), where k is Boltz- 
mann's constant, T « 290 K the ambient temperature, 
and where (•) denotes the average over many realizations 
of the Gaussian noise. The quantitative values of the re- 
maining model parameters 7, q* , and a are provided in 
the Appendices A.2-A.4. 

The model (J7J follows the standard approach for de- 
scribing electrokinetic effects and thus neglects non-linear 
terms of cubic and higher order in the external electri- 
cal field E(f,t), which would, in particular, take into 
account local changes of E(f, t) due to the presence of 
the particle. Likewise, hydrodynamic interactions are ne- 
glected. This approximation is justified in the present 
case of electrophoretically driven particle motion, be- 
cause in electrophoresis hydrodynamic interactions are 
essentially screened [40(. Remaining hydrodynamic ef- 
fects, in particular from interactions of the particle with 
the channel structure, are approximately taken into ac- 
count by determining the friction coefficient 7 from the 
particle diffusion in the microstructure, see Appendix 
A. 2. Since the friction forces are dominated by hydro- 
dynamic interactions with the bottom and top walls of 
the microfluidic channel (which are isotropic in x and y 
direction) , a possible orientation dependence of the fric- 
tion force in the x-y plane for particle movements parallel 
or perpendicular to the obstacle walls is expected to be 
negligibly small, and is thus omitted in our model. Ulti- 
mately, all these approximation are justified by the good 
agreement between the theoretical predictions and the 
experimental results (see, e.g., Figs. 5a and 6). 



C. Separation principle 

As already mentioned we want to address the prob- 
lem of separating an arbitrarily selected subset of particle 
species from a mixture of many different particle species. 
Obviously, the "optimal" realization of such a separation 
is to move the selected subset into one direction along the 
microstructured channel, while, at the same time, all the 
remaining species migrate into the opposite direction. In 



other words, the average particle velocity in x direction 



lim 

t— ¥00 



(<t) ~ *(0)) 
t 



(8) 



assumes different signs for the two sets of particle species. 
In the following, we devise the theoretical idea how this 
task can be accomplished with our microfluidic device 
from Fig. 1 by applying a periodic sequence of combined 
DC- and AC-voltage pulses. 

To keep things as simple as possible, we focus on 
cases where the DC- voltage mainly determines the elec- 
trophoresis, while the AC-voltage predominantly rules 
the dielectrophoresis. More precisely, we restrict our- 
selves to DC- and AC-voltages in ((TJ) such that practi- 
cally only the "systematic" DC-voltage C/dc contributes 
to the electrophoretic force q*E(r,t), while the effects 
of the "oscillating" AC-driving Uac sin(cji) occur on 
time- and length-scales, which are of no relevance for 
the effects we are interested in, and practically average 
out to zero during one oscillation. On the other hand, 
combining Eqs. (JTJ, ([2) and ([3]) we see that the time- 
averaged dielectrophoretic force aV£ 2 (r,t)/2 is propor- 
tional to aVE 2 (r)(U 2 c + U 2 AC )/U 2 = aVE 2 {r){\ + 
Ul c /Ul c )(U 2 AC /U 2 ). For |J7 DC | « Wkc\ and sufficiently 
large C/ac, the DC-contribution to the dielectrophoretic 
force aVE 2 (r,t)/2 can always be safely neglected, while 
the net effect of the AC-contribution over one oscilla- 
tion will still be significant. In other words, electrophore- 
sis is governed by the DC- and dielectrophoresis by the 
AC-voltage. In practice, these conditions are met by re- 
stricting ourselves to relatively small DC-voltages (a few 
volts in the experiment), while Uac may be substantially 
larger (several hundred volts). The AC frequencies lj/2it 
are restricted to the range of several hundreds up to thou- 
sands of Hz (a frequency of ui/2ir = 100 Hz turned out 
to be optimal for performing the experiments), thus also 
guaranteeing that the quasi-static approximation for a is 
valid (see also Sec. 2.1). 

From Fig. 2 we can infer that the electrophoretic forces 
alone would generate a net particle motion along the x- 
axis. On the other hand, the dielectrophoretic forces give 
rise to "particle traps" (attractors, deriving from local 
minima of the dielectrophoretic potential ([5])). Neglect- 
ing thermal noise effects for a moment, we are thus faced 
with a competition between directed motion (governed 
by the DC-voltage) and trapping (governed by the AC- 
voltage). Depending on the relative magnitudes of C/dc 
and C/aCj either the directed particle motion "wins" or 
the particle stays trapped. As a result, the entire pa- 
rameter space spanned by the DC- and the AC-voltage 
- henceforth denoted as "DC- AC-voltage plane" - can 
be divided into regions with steady particle motion and 
regions where the particle stands still, see Fig. 3a. Ob- 
serving that the electrophoretic forces are proportional 
to q*UuG an d the dielectrophoretic forces proportional 
to a U AC , the borderlines between those regions are gov- 
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FIG. 3: Theoretical concept, (a) Qualitative sketch of the 
parameter space spanned by the DC- and AC-voltage with n 
moving and TV — n trapped particle species for TV = 3 (see 
main text). The particle is trapped as long as the combined 
deterministic electrokinetic effects result in potential minima 
close to the obstacles (see Fig. 2b and main text), i.e. at the 
border between two different trapping regions (indicated by 
the solid lines), the potential minima just disappear for one 
of the particle species. Thermal noise effects are neglected 
in that picture; they would result in a finite trapping time 
which for that particle species approaches zero close to bor- 
der, but grows exponentially with the depth of the potential 
minima when moving away from the border into the trapping 
region towards larger Uac or smaller Uuc ■ For negative Uac , 
the regions are symmetric (see @) and hence omitted in the 
present plot. For negative Uuc, Eq. ([9]) implies an analogous 
symmetry apart from a different proportionality constant c. 
(b) Numerical solutions of © for particles with diameters 
1.1 /jm (black), 1.9 /jm (red), and 2.9 /jm (blue), moving in 
a periodically continued array of ratchet shaped obstacles as 
depicted in Fig. 2. The external driving is specified by Eqs. 
|T]), (J2J) and by Table 1, protocol a. Remaining parameter 
values and simulation details: see Appendices A.2-A.4 and 
A.5. 



erned by the equation 

Ul c = c\q*Uuc/a\ 



(9) 



with a proportionality constant c which depends on the 
details of the microstructure ( "forming" the spatial vari- 
ations of the electric field), the particle size, the sign of 
a, and - in case of asymmetric (ratchet shaped) obstacles 
- on the direction of the electrophoretic forces and thus 
on the sign of q*Unc, but otherwise is independent of q* , 
a, C/dc, and Uac- 

The main consequence is that for two particles species 
with different ratios q*/a or different sizes, the bor- 
derlines ([9]) between motion and trapping are different. 
(Here and in all what follows, we tacitly focus on the case 
that q* and a have the same sign for all species). As a 
consequence, the DC- AC-voltage plane now divides into 
regions with (i) both particle species moving, (ii) one 
moving and the other trapped, and (iii) both trapped. 
Likewise, for TV > 2 particle species, the DC- AC-voltage 
plane exhibits regions with n moving and TV — n trapped 
species for any n = 0,1, 2,..., TV. For an example with 
TV = 3, see Fig. 3a. 

By exploiting the existence of these regions, we will 
now construct a "voltage protocol" during which any 
particle of species k is on average displaced by Axk in 



x-direction, where k = 1, ..., TV labels the species and the 
Axfc are arbitrary (positive or negative) multiples of the 
lattice constant in x-direction. By repeating the corre- 
sponding protocol periodically in time, the average par- 
ticle velocity (jH)) of species k assumes the value Axk/r, 
where r denotes the duration of one period of the voltage 
protocol. By suitably choosing sign and modulus of the 
Axk, it is thus possible to make an arbitrarily selected 
subset of species move oppositely to all the others. In 
other words, any preset separation task can be accom- 
plished by properly tailoring the externally applied DC- 
and AC-voltage. 

Given the Axk for fc = 1, ..., TV, the announced proto- 
col is composed of TV + 1 time intervals of lengths At n , 
where n — 0,1, 2,..., TV and r = J2n=o^ n - Through- 
out the n-th time interval, the DC- and AC-voltages are 
kept at constant values U^q and U A( l, respectively. As 

seen above, the actual values of [/pg and U A q can and 
will be chosen so that n particle species are moving and 
TV — n are trapped. Furthermore, the indices k = 1, ..., TV 
of the particle species can be chosen without loss of gen- 
erality so that those with k < n are moving during time 
interval n and those with k > n stay trapped. Put dif- 
ferently, particle k (i.e. a particle of species k) is mov- 
ing during the time intervals with n > k and is trapped 
when n < k. Within the respective regions of the DC- 
AC-voltage plane (see above), UjQ and U A g can still be 

(n) 

freely chosen. In particular, U^q may still be of either 
sign, and likewise the transport direction of the moving 
species. 

We first focus on particle TV, which is only moving 
during the last time interval. Hence, Atjsr, U^ c , and 

U AC ', can (and must) be chosen so that particle TV cov- 
ers a distance Axn during this time interval TV. During 
the same time interval, all the remaining particles with 
k = 1,...,TV — 1 are also moving. Their respective dis- 
placements may be positive or negative and are denoted 
by Sxk- Once we made our choice of Ai^y, 17q C , and 

U AC ' also the traveling distances Sxk are uniquely fixed 
and can be calculated quantitatively. 

During all the remaining time intervals with n = 
0, ..,TV — 1, particle TV stays trapped and thus its to- 
tal displacement amounts to Axjy, as desired. On the 
other hand, the particles with k = 1, ..., TV — 1 now have 
to exhibit displacements by Ax' n = Ax n — Sx n during 
the remaining time intervals with n — 0, .., TV — 1. Along 
exactly the same line of reasoning, the existence of ap- 
propriate values of At n , U^l, and U A( l now follows by 
way of induction for n = N — 1,TV — 2, ..., 1. 

The left-over time interval with n = may appear su- 
perfluous since all particles are trapped. It is, however, 
essential for the separation procedure to work efficiently. 
In practice, there are various complications and perturba- 
tions of the above idealized scheme, most notably thermal 
noise and transient effects: after each "jump" of E/dc and 
L/aCi every particle needs a certain "relaxation" time to 



adapt its motion to the new circumstances (approach the 
new attractor of the dynamics) , especially when switch- 
ing from "trapped" to "moving" or vice versa. All these 
effects would prevent a reliable time-periodic particle mo- 
tion upon periodically repeating the protocol without the 
"all-trapping" time interval n = 0. Moreover, it guaran- 
tees that the particles are indeed trapped before species 
by species is released from the traps during the next cycle 
of the protocol. 

We emphasize that the only indispensable prerequisite 
of the above procedure is the existence of the different 
regions in the DC-AC-voltage plane (Fig. 3a), which in 
turn represents the generic case. Different electrophoretic 
mobilities of the particle species or any other special sys- 
tem or particle properties are not required. Accordingly, 

(n) (n) 

the parameters At n , C/f) C , and C/^ c are not yet uniquely 
fixed and thus can still be optimized with respect to the 
experimental constraints and uncertainties. Along these 
lines one typically ends up with alternating signs of U^q 

and decreasing values of \U^\ (the sign is irrelevant). 
Moreover, in some cases one finds that less than N + 1 
time-intervals may already do the job, see Fig. 3b and 
Sec. 3. 

In principle, by the above strategy one may separate 
particle species of arbitrary sizes or q* jot ratios in a 
rather broad dynamical range. In particular, increas- 
ing noise effects for smaller particles can be kept under 
control by applying larger voltages and hence also en- 
hancing the deterministic effects. Likewise, reasonably 
small particle-to-particle variations within one species es- 
sentially lead to a shrinking of the pertinent regions in 
the DC-AC-voltage plane (cf. Fig. 3a) and hence are still 
admissible. However, too large variations lead to a dis- 
appearance of the regions and thus to a breakdown of 
our approach. Analogous considerations apply to devi- 
ations from a strictly periodic microstructure. On the 
other hand, a mixture of particles with continuously dis- 
tributed sizes or q* jot ratios may still be sorted into cer- 
tain subclasses, but a truly continuous separation is not 
possible. 

Our general scheme would in principle also work for 
symmetric obstacles. Ratchet shaped obstacles have the 
advantage of different c-values for positive and negative 
DC-voltages in §§§, admitting much more efficient and 
robust separation in practice |27| . Likewise, the above 
mentioned transients after each "jump" of Uuc and Uac 
can be exploited to further improve the efficiency of the 
procedure, especially when taking into account thermal 
noise effects. We finally note that somewhat similar con- 
cepts in the context of ratchet effects and particle sorting 
have been utilized previously in Refs. [4l| - |43l |. 



III. RESULTS 

To quantitatively verify these heuristic ideas, we per- 
form detailed numerical simulations of the model (|7|). see 



Appendix A. 5. As a first example, we consider N = 3 
particle species. Fig. 3a illustrates the corresponding re- 
gions in the DC-AC-voltage plane with n — 0, 1, 2, 3 mov- 
ing and N — n trapped species. Fig. 3b exemplifies the 
separation of one species from all the others: During an 
initial time-interval of 10s (see Table 1, protocol a), all 
particles are trapped. During the second interval of 50 s, 
two species are moving forward and one remains trapped. 
During the last interval of 12 s, all move, but now into 
the negative direction. Upon periodically repeating the 
protocol (two periods are shown), two species acquire a 
positive and one a negative average velocity. (An inter- 
val with only one moving species is not needed for this 
specific separation task.) By suitably tailored protocols 
(see Table 1), also any other of the three particle species 
can be made to migrate oppositely to the others, as will 
be experimentally demonstrated below. 

Before doing so, a second illustration with N = 5 par- 
ticle species is presented with Fig. 4. Proceeding as de- 
tailed above, we have tailored voltage protocols for three 
separation tasks which appeared to us particularly in- 
teresting and non-trivial: Fig. 4a demonstrates the sep- 
aration of the middle sized particles from the two big- 
ger and the two smaller species. Fig. 4b shows how the 
largest and smallest particle species can be made to move 
oppositely to the three middle sized species. The most 
challenging case is shown in Fig. 4c, namely the disen- 
tanglement of the particles with 1 /im, 2 //m, and 3 /jm 
diameter from those with 1.5 ^m and 2.5 /um diameter. 

Closer inspection of Figs. 4a-c reveals that the particle 
motion does not always precisely follow the displacement 
pattern anticipated in the construction of the voltage pro- 
tocol as described above (Sec. 2.3), but rather shows sig- 
nificant variations between different periods of the pro- 
tocol. For instance, the 2 /im (green) particle in Fig. 4c 
moves in four (out of the six) voltage steps during the first 
period and only in three during the second period. Like- 
wise, the 2.5 ^tm (blue) particle shows a net displacement 
of about 30 /xm during the first cycle but performs no 
net movement in the second cycle. These variations are 
due to thermal noise effects, which have a strong impact 
on the escape process out of the dielectrophoretic traps, 
because these are located close to stagnation points of 
the DC electric field (see Fig. 2). Due to this variabil- 
ity of the net displacement per voltage cycle, particles 
of the same species are spread out along their direction 
of systematic motion farther than by pure thermal dif- 
fusion. Such enhanced dispersion notwithstanding, effi- 
cient particle fractionation based on transport in oppo- 
site directions is demonstrated by Fig. 4d for all three 
separation tasks from Figs. 4a-c. Note that even within 
the subsets of particle species, that have been selected to 
migrate in one specific direction, different migration ve- 
locities for different species and thus (partial) separation 
may be observed (e.g. in the first panel of Fig. 4d). This 
"side-effect" results from performing the separation by 
less than N + 1 steps in the voltage protocol, so that the 
displacements A.Xk have not been determined and opti- 



TABLE I: Voltage protocols for the simulations in Figs. 3b, 6 and the experiments in Fig. 5. Listed are the voltages Uuc and 
Uac in volts and the time-intervals At n in seconds during which they were kept at those constant values. The entire protocol 
U(t) in {TJ is obtained by periodically cycling through the listed sequence of time-intervals. In all cases, the AC-driving 
frequency in fT]) is set to uj/2-k — 100 Hz. 



Protocol a 


Protocol b 


Protocol c 


Atn 


Uuc 


Uac 


Atn 


Udc 


Uac 


Atn 


Uuc 


Uac 


10 





400 


10 





400 


10 





400 


50 


-15 


320 


90 


+ 15 


400 


40 


+15 


360 


12 


+10 





25 
10 


-16 
+ 10 


330 



12 


-10 






(a) 




(c) 




u 500 

r-, 20 



I ' I 




20 30 40 50 

»[s] 




(b) 



-100 

—.1000 

500 
a 

> 20 

a 
a 

& -20 

I 

(d) 0.4 

£ 0.3 

1 0.2 
I 0.1 





40 60 

<[s] 




I ' I 




10 20 30 40 50 60 
t [S] 





-2000 


2000 


4000 


& 0.4 




| ' 


' 


' : 


1 0.3 
1 0.2 
1 0.1 


L J 


l 


A 


J 


r 



S- 0.3 

1 0- 2 

2 0.1 

a. 






-2000 2000 

x-position Qim] 



FIG. 4: Numerical solutions for 5 particles with diameters 
lpm (black), 1.5 pm (red), 2pm (green), 2.5 pm (blue), and 
3pm (magenta), (a): Separation of the 2pm particles (mid- 
dle size) from all the other, (b): Separation of the smallest 
and largest (lpm and 3 pm) particles from all the others (3 
intermediate diameters) . (c) : Separation of the particles with 
1 /im, 2 pm, and 3 pm diameter from those with 1.5 pm and 
2.5 pm diameter. The lower panels specify the voltage proto- 
cols used to achieve the separation tasks shown in the respec- 
tive upper panels. In each case, two periods of the voltage 
protocols are shown. The explicit algorithm to devise those 
protocols is given in the main text. Further details and pa- 
rameter values are as in Fig. 3b, see also Appendices A.2-A.5. 
(d): Spatial distributions of 100 particles of each species ob- 
tained after simulating 100 periods of the voltage protocols. 
The upper panel corresponds to the separation task from (a), 
the middle panel to (b), and the lower panel to (c). 



mized independently for all the species in the subset (see 
the algorithm for constructing the voltage protocols, as 
described in Sec. 2.3), and from the fact that the devia- 
tions from the anticipated displacement patterns depend 
on particle properties and arc thus different for different 
particle species within one subset. 

Next, we turn to the experimental test of our theoret- 
ical concepts by means of the setup as detailed above in 
Fig. 1 and in Appendix A.l. Fig. 5 summarizes our main 
findings, demonstrating that for a colloidal suspension 
containing three different species of polystyrene beads it 
is indeed possible to control the transport directions to 
the extent that any one of the three species can be made 
to propagate oppositely to the other two species. 

Finally, Fig. 6 shows numerical simulations, obtained 
along the very same lines as in Fig. 3b, except that the 
so generated numerical data have now been further pro- 
cessed exactly as the experimental data in Fig. 5a. These 
numerical results in Fig. 6 compare very well with the 
experimental counterparts from Fig. 5a. The fact that 
the velocity spreads in Fig. 5a are not significantly larger 
than in Fig. 6 implies that these experimental spreads 
are mainly caused by the thermal noise effects discussed 
above and not by variations of q* , a, and 7 for beads of 
the same species. The minor remaining differences be- 
tween experiment and simulations are probably due to 
nonlinear (induced-charge) electrokinetic effects [4J, |45| , 
which we have neglected in our model (0. Furthermore, 
it can not be fully excluded that (small) Joule heating ef- 
fects [46| affect the particle motion, although we do not 
observe any of the typical indicators for Joule heating 
during the experiments, like vortex-like fl ow p atterns at 
the constrictions created by the obstacles [47], |48| or the 
generation of gas bubbles. 



IV. CONCLUSIONS 

In summary, we have put forward a new concept for 
particle sorting in a microstructured ratchet device by ex- 
ploiting a subtle competition between electrophoretic and 
dielectrophoretic forces. Within one and the same device, 
by tailoring the applied voltage protocol, one or several 
arbitrarily selected particle species can be made to prop- 



agate oppositely to all the other particle species. The 
choice of those particular species can be easily changed 
and adapted at any time, even while the experiment is 
running, by simply adapting the applied voltage protocol. 

Such versatile particle transport is achieved by exploit- 
ing all "elementary ways of symmetry breaking" in the 
context of ratchet effects [Irtl , namely spatial symmetry 
breaking, temporal (or dynamical) symmetry breaking, 
and application of a bias force. Any demand to maintain 
certain symmetries would amount to an unnecessary con- 
straint, whereas their combined breaking allows for the 
required flexibility and efficiency |27| . 

Our new separation concept has been devised theo- 
retically and realized experimentally within a particle 
density regime, in which particle interactions other than 
excluded volume effects are negligible and in which the 
number of dielectrophoretic particle traps is larger or 
comparable to the total number of particles in the mi- 
crofluidic device. 

In the future, experiments with biological complexes 
ranging from single biopolymers, like DNA or proteins, 
to whole cells, as well as non-biological objects, like 



£> 0.4 


» 


i i 


i i i i i 


exp. data " 


J 0.2 

o 

& o- 


n_ 


4-0- 


, -rvft 


v^v 



-1.0 -0.5 0.5 1.0 1.5 2.0 2.5 



& ° 4 

I 0.2 

o 

S-H 

& o 



(b) 



1 — i — i — ' — i — ' — r 



exp. data 




-1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 





(c) 


exp. data - 


£> 0.2 










n r n - 


S 0.1 


■ 


L 










s 

& n- 


i-i|i~i | ,0-, .r— ^ ■ — ■ i 


h r-ri n ft. 


M 



-2.0 -1.5 -1.0 -0.5 0.5 1.0 1.5 2.0 
velocity [iim/s] 

FIG. 5: Experimentally measured velocity histograms for 
N — 3 particle species (polystyrene beads with diameters 
1.1 (jm (black), 1.9/tm (red), and 2.9 /xm (blue)), simulta- 
neously moving in the structured microfluidic ratchet device 
from Fig. 1. (a)-(c) show the results for the three different 
time-periodic voltage protocols from Table 1. For each his- 
togram, the displacement of 100 beads in x-direction during 
one time-period r of the driving voltage has been measured 
and interpreted as their net velocity during that driving pe- 
riod, i.e. an average over the histogram for a specific particle 
species gives the average velocity ([5| for that species. In (a), 
the 2.9 /mi beads move oppositely to the other two species 
with very high reliability, and analogously for the 1.9 /mi 
beads in (b) and the 1.1 /mi beads in (c). 
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experiment from Fig. 5a according to {TJ, (JTJl, and Table 1. 



carbon-nanotubes, can be envisioned. In spite of the 
fact that dielectrophoretic forces are theoretically ex- 
pected to strongly decrease with particle size (cf. Eq. 
©), the feasibility of manipulating single proteins by 
electrodeless dielectrophoresis has already been experi- 
mentally demonstrated e.g. in |49l - [51| . Moreover, sepa- 
ration may be performed not only according to size but 
also according to molecule conformation or particle shape 
[52j and thus represents a considerable advance beyond 
conventional separation techniques, like (pulsed-field) gel 
electrophoresis. A further interesting direction for future 
research are hydrodynamic and other interaction effects 
as well as possible modifications of the separation mech- 
anism, arising for substantially higher particle densities 
than in our present experiments. 

This work was supported by Deutsche Forschungsge- 
meinschaft under SFB 613 and RE1344/5-1 



Appendix A: Materials and Methods 
1. Experimental details 

The entire microdevice (see third paragraph and 
Fig. 1) was cast with poly(dimethylsiloxane) (PDMS), 
mounted on a PDMS coated coverslip after a treatment 
with an oxygen plasma for 30s, and coated with a com- 
mercial triblock copolymer (F108, 500 /iM) for 2 hours 
53]. Spherical polystyrene beads (Molecular Probes) 
with diameters 1.1 /mi, 1.9 /mi, and 2.9 /mi at a concen- 
tration of 50 - 500 particles per species and nanoliter were 
immersed in a 10 mM phosphate buffer (pH 8.3) con- 
taining 200/xM Tween 20. Hence, observable hydrody- 
namic or other interaction effects were always negligibly 
rare. For a more detailed description of our PDMS chip 
fabrication, par ticle and microstructure pre-treatment, 
we refer to [13, H]|. Note, that the F108 coating of 
the microstructure strongly suppresses electroosmotic ef- 
fects [53| , as well as adhesion of the polystyrene beads to 
the PDMS surfaces (in the experiments less than 1% of 
the particles got stuck to the PDMS structure). Parti- 
cle tracking by video microscopy was carried out with 
a standard fluorescence microscope. We recorded the 
motion of the particles within a field of view of about 
825 /mi x 310 /tm with a frame-rate of 5 fps. U sing the 
software Image J 1.43 (with the plugin MtrackJ) [56ll57|. 
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particle trajectories were extracted from the recorded 
videos with a spatial resolution of 0.6 /im (both in x and 
y direction) and a temporal precision of 0.2 s. The three 
different particle species were easily distinguishable by 
their size and/or fluorescent labeling (Fig. lc). The his- 
togram plots in Fig. 5 were obtained from the trajecto- 
ries by determining the particle displacements during one 
period of the voltage protocol; the so obtained velocities 
have a relative error below 10%. 



2. Friction coefficient 

For a spherical particle with diameter d in a fluid 
of viscosity r\ and infinitely far boundaries, the fric- 
tion coefficient 7 in Eq. (J7]) could be approximated by 
Stokes' formula ^stokes = 67rry(d/2). Here, we go one 
step beyond this approximation, accounting very roughly 
for the hydrodynamic interactions with the microstruc- 
ture by experimentally determining the diffusion coef- 
ficient D and then inferring 7 from Einstein's relation 
D = fcT/7. Those experimentally estimated values 
were D = 0.202 ± 0.002 /im 2 /s for the 1.1 /xm particles, 
D = 0.131 ± 0.005 /im 2 /s for the 1.9 /im particles, and 
D = 0.082 ± 0.007 /im 2 /s for the 2.9 /im particles. The 
resulting 7- values are close to 1.7 • 7 stokes in all three 
cases. For the five particle species simulated in Fig. 4 
we thus assume along the same line of reasoning that 

7 = 1.7 • JStokes- 



3. Effective charge 

The effective charge q* from (J3J) depends in a very com- 
plicated way on the particle material, the buffer solution, 
and the surface coating of particle and microstructure 
|33l l3a . |58| . In principle, one would expect that parti- 
cles which only differ by their diameter have the same 
zeta potential, and hence the effective charge @ would 
inherit the linear dependence on the particle diameter 
from the friction coefficient 7. In practice, even nomi- 
nally identical particles but purchased at different times 
or from different companies are found to exhibit differ- 
ent electrophoretic mobilities due to somewhat differing 
production conditions and also due to ageing effects of 
the surface coating (induced by hardly controllable (bio- 
chemical reactions at the particle- fluid interface). In 
view of these practical problems, we have directly es- 
timated q* by fitting the model ([7]) to the experimen- 
tally observed particle motion, resulting in q* = —0.055 • 
10~ 15 C for the 1.1 /tin particles, q* = -0.17 • 10~ 15 C for 
the 1.9 /im particles, and q* = -0.50 • 10~ 15 C for the 
2.9 /im particles. 

For the five particle species simulated in Fig. 4, we as- 
sumed - in the absence of experimental data - the above 
mentioned idealized theoretical situation of identical zeta 
potentials and thus identical Helmholtz-Smoluchowski 
velocities Q. In the concomitant linear dependence of 



q* on the particle diameter, the proportionality constant 
was chosen so as to yield a net velocity of 6 /im/s in 
x-direction when driven by a purely static voltage of 
Ubc = 10 V. This velocity is close to those we observed 
for the three experimental particle species. 

We remark that the electroosmotic fluid flow v{r) of 
the ionic liquid exhibits a remarkable "similitude" to 
the electrophoretic Helmholtz-Smoluchowski relation ([3]), 
namely v(r) = —E(r)e(' /rj. where £' is the zeta potential 
of the microstructure, see [58|, [59j and further references 
therein. The effect of this electroosmotic flow on the 
particle can be readily incorporated into the above men- 
tioned effective electrostatic force q*E(r) by replacing £ 
in the definition Q by £ — £'. In our actual experiment, 
those electroosmotic effects are, however, strongly sup- 
pressed by the F108-coating of the microstructure (see 
above) . 



4. Effective polarizability 

While the particle polarizability is in general frequency 
dependent, experimentally we always work in the quasi- 
static regime (uj/2tt = 100 Hz), and hence our a always 
refers to the zero frequency limit. Invoking the basic for- 
mula © and approximating the medium's permittivity 
e by that of water, e = 80eo (eo being the vacuum dielec- 
tric constant), one finds that a = —7.4 • 10 -28 Fm 2 for 
the 1.1 /im particles, a = —38 • 10~ 28 Fm 2 for the 1.9 /xm 
particles, and a — — 136- 10~ 28 Fm 2 for the 2.9 /im parti- 
cles. The negative sign of a is in agreement with previous 
experiments under comparable conditions and has been 
tacitly anticipated in Fig. 2b. 

However, the basic formula © neglects, among others, 
the polarizability of the counter-ion cloud (Debye layer) 
around the dielectric sphere, which, in turn, is much more 
difficult to estimate theoretically. Like for the above dis- 
cussed effective charges, we therefore estimated the rel- 
evant effective polarizabilities by directly fitting the ex- 
perimentally observed particle motion to the theoretical 
dynamics from Eq. (JTJ), resulting in a — —4.0 • 10~ 28 
Fm 2 for the 1.1 /im particles, a = —27 • 10 -28 Fm 2 for 
the 1.9 /im particles, and a — —126 • 10~ 28 Fm 2 for the 
2.9 /im particles. The slight systematic deviations from 
([6"J) are remarkably small. For this reason and due to the 
lack of experimental data, the polarizabilities for the five 
particle species simulated in Fig. 4 were chosen according 
to formula ^ without any further corrections. 



5. Simulation details 

For numerical simulations, the model (O was dis- 
cretized according to the Euler algorithm with a time 
step of dt = 10 /is. The electric field E\ (r) and the di- 
electrophoretic "field" V-E 2 if) was computed on a spatial 
grid with step sizes 0.02 /xm within one basic "unit cell" of 



the microstructure around an obstacle (cf. Figs. 2a, b) by 
using the COMSOL Multiphysics software suite. This 
was done once before performing the dynamic simula- 
tions based on ([7]). The actual electrophoretic and dielec- 
trophorctic forces at the current particle position during 
the simulations were determined by periodically contin- 
uing the fields within this "unit cell" (in x and y di- 
rection), by interpolating linearly between the relevant 
grid points, and by taking into account ((2|). The term 
F(f) in ([TJ, representing idealized hard- wall interactions 
with the microstructure, was treated according to the so- 
called rejection scheme [6(| ■ The triangular posts of the 



microstructure were approximated by the shape shown 
in Fig. 2. 

To obtain the simulation results of Fig. 4, 100 indepen- 
dent realizations were computed for each particle species 
over 100 periods of the respective voltage protocol (the 
Figs. 4a-c show one typical realization for each particle 
species during the first two driving periods). In Fig. 6, 
the results of 100 independent realizations for each par- 
ticle species are shown, evaluated over one period of the 
voltage protocol (in accordance with the processing of 
the experimental data). 
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